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FOREWORD 


This  report  was  prepared  by  the  Hypersonic  Research  Laboratory  of 
the  Aerospace  Research  Laboratories,  Air  Force  Systems  Command, 

United  States  Air  Force,  under  Project  7064  entitled,  "High  Velocity  Fluid 
Mechanics.  " 

The  results  given  in  this  report  are  a  by-product  of  a  study  of  the  effects 
of  Mach  number  gradients  prevalent  in  wedge-shaped  or  conical  nozzles 
employed  for  hypersonic  wind  tunnel  testing. 


ABSTRACT 


The  effect  of  high  Mach  number  flows  on  test  bodies  can  be  estimated 
reasonably  well  by  using  the  Newtonian  approximation.  If  the  flow  is  like  a 
source  flow  then  the  effects  of  the  ambient  flow  gradients  on  the  flow  para¬ 
meters  in  the  shock  layer  can  be  obtained  by  modifying  the  Newtonian  formula.: 
Freeman's  theory  is  adopted  in  the  present  study  and  changes  in  surface 
pressure  and  shock  detachment  distance  of  a  two-dimensional  blunt  Lody  are 
obtained  by  using  a  perturbation  scheme  Based  on  the  ratio  of  the  axial 
distance  of  a  station  on  the  body  from  the  nose  (which  is  the  origin  of  the 
coordinate  frame)  to  the  distance  of  the  "  source"  location  from  the  origin 
as  a  perturbation  parameter,  first  order  corrections  to  the  surface  pressure 


and  the  shock  detachmert  distance  are  derived. 


INTRODUCTION 


The  plane  and  axisymmetric  hypersonic  flow  past  blunt  bodies  have 
been  investigated  as  an  inverse  problem  (shock  shape  being  given)  by 
Chester^  and  Freeman^.  The  fluid  is  assumed  to  be  in  thermodynamic 
equilibrium  locally,  and  viscosity  and  heat  conduction  are  neglected.  The 
governing  flow  equations  for  a  perfect  gas  are  expanded  in  powers  of 
(y -  l)/(y  +  1)  and  Mi2,  where  y  is  the  ratio  of  the  specific  neat?  and  Mq, 
is  the  Mach  number  of  the  free  stream.  The  first  approximation  correspond¬ 
ing  to  y  -  1  and  Ma>=  ®  yields  a  pressure  formula  identical  with  a  formula 
(3) 

given  by  Busemann  .  Busemann  pressure  formula  is  often  referred  to  as 
"  Newton- Busemann  pressure"  or  "  Newtonian-centrifugal  pressure.  "  In 
regions  of  the  flow  where  pressure  is  not  a  smaU  fraction  of  its  value  at  the 
nose,  the  theory  of  "thin  shock  layers"  of  Chester  and  Freeman  yields  con¬ 
sistent  expressions.  However,  at  some  distance  downstream  of  the  stafcnarion 
point,  Freeman's  theory  becomes  invalid  as  the  pressure  becomes  ?.  small 
fraction  of  its  value  at  the  nose.  A  "  singularity"  appears  which,  f example, 

on  a  sphere  will  be  located  at  a  central  half  angle  of  60°.  Schneidex'3 

(4) 

theory  removes  this  difficulty,  ana  ,  utvides  a  more  uniformly  valid  solution. 
The  thin- shock  layer  assumption  is  not  made  in  Ref.  4,  and  changes  in  velo¬ 
city  along  streamlines  is  permitted.  Schneider's  theory,  thereby,  improves 

Freeman's  analysis.  (Freeman  has  also  modified  his  theory  to  give  a 

(5) 

uniformly  valid  solution  in  the  region  of  the  singularity  .)  Freeman's 


i 


theory,  however,  describes  the  flow  conditions  near  the  stagnation  region 
reasonably  weii.  This  theory  is  applied  to  give  the  shock  detachment  distance 
when  the  flow  is  a  two-dimensional  source-type  flow  on  a  two-dimensional 
blunt  body.  This  condition  corresponds  to  the  flow  in  a  wedge  shaped  nozzle 
in  a  hypersonic  wind  tunnel.  The  shock  stand-off  distance  increases  from 
its  uniform-flow  value  due  to  Mach  number  gradient  in  the  wedge  shaped 
nozzle.  The  decrement  in  the  surface  pressure  due  to  the  gradient  effect  is 
also  obtained  by  following  a  perturbation  scheme. 
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ANALYSIS 


The  flow  in  a  wedge  shaped  nozzle  can  be  considered  to  be  a 
source-like  flow.  The  source  is  assumed  to  be  located  near  the  throat  of 
the  nozzle.  Any  test  body  in  the  wind  tunnel  exposed  to  a  source-like  I  yper- 
sonic  flow  will  experience  the  effects  of  the  Mach  number  gradients  which 
arise  due  to  the  angularity  of  the  flow.  Freeman's  theory  is  extended  to 
obtain  the  corrections  due  to  the  flow  angularities.  Analysis  of  the  adaption 
of  Freeman's  results  is  described  below. 

The  ideal  source  flow  in  two- dimensions  leads  to  the  equation 

Pl  qj  r  =  const, 

where  pi  and  qt  are  ambient  density  and  velocity  respectively.  The  term 
r  represents  the  radial  distance  of  the  reference  station.  If  the  ambient 
values  are  assumed  to  be  quantities  perturbed  about  their  respective  values 
(indicated  by  barred  quantities)  at  the  nose,  we  write  the  equation  as 

(pi  +  Api)  (qi  +  AqJ  (rQ  +  x')  =  piqir0  . 

Here  r  is  approximated  by  r0  +  x^  where  r0  represents  the  distance  of 
the  nose  from  the  throat,  and  x f  is  the  axial  coordinate  measured  (from 
the  nose)  of  a  reference  station  on  the  body.  The  approximation  rQ  +  x*ss  r 
is  reasonable  as  we  are  concerned  with  a  region  in  the  vicinity  of  the 
stagnation  point.  Equation  (2)  can  be  written  as 


it  being  assumed,  that  higher  powers  of  x*/  rQ  are  small  quantities  and  are, 
therefore,  negligible.  Since  in  high  Mach  number  flows  the  relative  change 
in  axial  velocity,  Aqx  /qj  ,  is  small  compared  to  the  relative  change  in 
density,  &f>  /px,  Eq.  (3)  becomes 

I21  =  -  -  •  (4) 

Pi  ro 

Also  from  the  relation 

Pi  ~  ypi 

we  write 

^  .  (5) 

Equations  (4)  and  (5)  are  now  adopted  in  deriving  the  corrections  due  to  the 
ambient  flow  nonuniformi  ies. 

The  Newtonian- centrifugal  expression  for  the  pressure  behind  the 
shock  is  given  by  (Ref.  2) 

~  ®in2  a  (x)  +  j  ^  (t)  cos  cr  (t)  dt  .  (6) 

x 


The  coordinate  system  is  body  oriented,  x  denoting  the  curvilinear  coor¬ 
dinate  along  the  body  surface  and  y  denoting  the  normal  coordinate  to  the 
body.  The  terms  in  Eq.  (6)  represent  the  following  quantities: 

Pi  =  ambient  density  in  the  nonuniform  flow 
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pressure  in  the  shock  layer  at  the  pcint  of 
intersection  of  the  streamline  0  -  const,  and 
the  normal  (to  the  body  surface)  at  x.  The 
reference  streamline  0S  intersects  the  shock 
at  a  point  Ps  whose  body-coordinate  is  x. 
Pressure  at  Pg  behind  the  shock  is  given  by 
the  Rankine-Hugoniot  condition 

ambient  velocity  in  the  nonuniform  case 

angle  which  the  ambient-flow- streamline  in  the 
nonuniform  case  makes  with  the  shock  at  Pg 

angle  of  inclination  of  the  ambient  streamline 
in  the  nonuniform  flow  case  with  the  axis  of 
symmetry 

geometric  scale  factor  which  is  equal  to  unity 
in  the  two-dimensional  case 

the  curvature  of  the  shock  (and  hence  of  the 
body  according  to  the  thin  shock  layer  assump¬ 
tion. )  This  is  also  given  by  -  <f/ (x)  where 
0  (x)  refers  to  the  body  slope  (approximately 
equal  to  the  slope  of  the  shock)  at  the  point  x. 


Following  Freeman's  notation,  we  introduce  a  quantity  7}(x)  which  refers 
to  the  distance  of  a  point  on  the  body  from  the  line  of  symmetry.  Therefore, 
♦  (t)  as  7j(t)  in  the  two-dimensional  case.  In  the  equations  hereafter,  we 
use  the  subscript  "  s"  to  denote  quantities  behind  the  shock. 

Introducing  the  perturbation  quantities  in  Eq.  (6),  we  get 


plg(x,{)  +  Api8(x,$)  .  i.  .  f*  . 

i — aoTT-V/  "au1V  =  sin  W- ^  +  P(x)  /  ^ 


(t)cos  (0  -  ji)dt  , 
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Pl  U1  [  PlSj  L  Pl  J  L  ’  2~ui  J  “  [8in^^  "  M(x)  COB  0fx)J 

+  f'tx)/' T7'(t)  [cos^(t)  +  M(t)  sin*  (t)]  dt.  (7) 

In  Eq.  (7),  it  is  assumed  that  angle  p  is  small  and  ^  /rQ.  Also 
powers  higher  than  unity  of  the  relative  changes  of  the  flow  quantities  are 
neglected  as  are  the  cross  products  of  the  quantities  Apj8,  APi  and  AUX  . 
Equation  (7)  can  then  be  written  as 

*  If1  TTZf  +  <l',x)  J  ’’<(t)  t‘(t)  8in * <*>  dt 

(8) 

-  p(x)  sin  2p  (x)  . 

We  note  that  ( x *  being  the  coordinate  of  the  reference 

point  on  the  body  measured  aljng  the  axis  of  symmetry  from  the  nose.  )  Also 

p(x)  =  V  (x)/ro. 

In  Eq.  (8),  the  term  pls  is  given  by  the  expression 

^L8  %’}—  ^  sin2  0  (x)  +  /(x)  ^  7?'(t)  cos  0  (t)  dt  . 

Pl 

It  can  be  seen  that  Eq.  (8)  gives  in  the  vicinity  of  the  blunt  stagnation  region  a 

surface  pressure  correction  which  is  approximately  one  half  of  the  pressure 

perturbation  for  a  sharp  edged  body.  This  is  in  qualitative  agreement  with  the 

(6) 

work  of  Gordon  Hall 

The  equation  for  the  streamlines  in  the  nonuniform  flow  case  is  given  by 


y  =  € 


/ - 

o  |sinz(7(x)  + 


.  sin  4>  (t)  sinzq(t)  dt _  ^ 

sin^(t)  cos  Oft)  dtj  |cos*a(t)  -  2elogj^^7^  |  z 


where  €  =  (y  -  1  )/(y  +  1)  (Y  being  the  adiabatic  exponent).  Further 
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*  sin2a(x)  +  0*  sin  0(t)  cos  (7  (t)  dt 


We  note  here  that  in  the  two-dimensional  case 
♦V)  ~  7)fy)  «  sin  0  (t) 

After  introducing  substitutions  and  approximations  as  were  done  earlier, 
we  observe  that  Eq.  (10)  becomes 

=  sinz0  (x)  -  fi(x)  sin  2  0(x) 

Pi  ui 


+  0*(x)  y  sin 


0  (s)  cos  0  (s)  ds  +  I  jji(s)  sin20(s)ds|  (11) 


We  note  further  that  in  Eq.  (9)  the  expression 


R1  (x,  t)  *  sin2  a(x)  +  0*  J  sin  0(s)  cos  a(s)  ds 


can  be  written  as 


R^x.t)  a  sin2  0(x)  -  H(x)  sin  2  0(x) 


+  0#(x) 


(s)  cos  0  (s)  ds  +  I  ji(s)  sin20(s)  ds 


In  view  of  Eqs.  (12)  and  (13),  Eq.  (9)  can  be  written  as 

_ fsin^  0  (t)  -  u(t)  sin  p  (t)  sin  2  0  (t)  1  dt 

v  R*  (x,  t)  [cos20  (t)  +  (1  sin  2<p  (t)  -  2  e  log  R1  (x  o)J  ^ 

The  equation  for  the  shock  is  obtained  by  replacing  £  by  x.  We  are 
now  able  to  evaluate  the  shock  "  stand-off"  distance  in  the  two-dimensional 
case.  We  note  that  for  small  values  of  t,  T)(t)~t,  and  0  »  U  /  2.  -  t/a, 
where  a  is  the  radius  of  curvature  of  the  nose.  With  these  approximations, 


7 


Eq.  (14)  becomes 


fcos3  —  cos—  sin  — 1  dt 

_ L _ .f° _ 5 _ _ _ 

R1  (x,  t)  [si“2“  +  7-  sin~  -  2  c  log  Rl  (x,  o)  2 


[15) 


where  fi(x)  ~  r]x/  rQ  ~  x/rQ  when  x  is  small, 
and 

R1  (x,t) 


2  X  1  2  t  x  .  2x  1  .  2  ^2. 

—  -  —  cos  —  -  sin  —  -  — 1  (x  -  t  ) 


3 

-  cos 
Z  a  Z 


a  r„ 


4ar 


8r 


f2x  . 

I —  si 

La 


2x  2x 

sin  —  +  cos  — 
a  a 


2t  2t  2t  I 

-  —  sin  —  -  sin  — 
a  a  a  J 


After  some  algebra,  Eq.  (15)  can  be  written  as 


y  =  € 


/*  ( cos3  — )  dt 

I  V  a/  _€_  I  F(x,  t)  dt 

bf  D  ro  D 


(16) 


where 


and 


“■(i 


2Xl  j 

COS  —  -  —  COS 

a  Z 


t\  /  2  t 

-  I  (sin  - 

a/\  a 


+  2  e  log 


( 


3  cos2 


2  \Y‘ 


2  cos 


(3  cos2  —  -  cos 
\  a  a 

t  sin—  os4  — 
a  a 


x2  ,  5x  .  2x  ,  a 

*7\  T~  +  "7*  sin  —  +  - 

t\  4a  4  a  8 


sin2  —  +  2  elog 


-  2  t  sin  —  cos2  — 
a  a 


2x  t2  t  2t 

cos  —  -  —  -  ~  sin  — 

a  4a  4  a 

a  .  2t  1 

'  I8m  T 

t  2 1  J 


r3  cos 


2  x 


The  first  integral  on  the  right  hand  side  of  Eq.  (16)  can  easily  be  evaluated, 
and  in  the  limit  as  x->0,  the  value  of  the  integral  (which  is  ae  /2  log(4/  3e)) 
gives  the  shock  stand-off  distance  in  the  uniform-flow  case..  The  second 
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integral  on  the  right  hand  side  of  Eq.  (16)  is  the  correction  to  the  stand-off 
distance  due  to  the  ambient  flow  nonuniformities.  The  integral  can  be 
evaluated  in  a  straight  forward  manner,  and  the  value,  in  the  limit  as  x->0; 
becomes 


*'  *  icj-  £1°8^  • 

We  can,  therefore,  write  the  shock  stand-off  distance  in  the  nonuniform  case 
as 


where 


and 


(17) 


(in  the  uniform  flow  case) 


(18) 


A8 

a 


,  (correction  due  to  the  ambient 
flow  nonuniformities) 


(19) 


Equation  (19)  gives  the  increment  in  the  shock  stand-off  distance  due  to 
angular  divergence  of  the  ambient  flow.;  The  increment  can  be  quite  significant 
when  the  radius  of  curvature  at  the  nose  of  the  test  body  is  large. 
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